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Abstract 

We suggest a generalization of the Lie algebraic approach for con- 
structing quasi-exactly solvable one-dimensional Schrodinger equati- 
ons which is due to Shifman and Turbiner in order to include into 
consideration matrix models. This generalization is based on repre- 
sentations of Lie algebras by first-order matrix differential operators. 
We have classified inequivalent representations of the Lie algebras of 
the dimension up to three by first-order matrix differential operators 
in one variable. Next we describe invariant finite-dimensional sub- 
spaces of the representation spaces of the one-, two-dimensional Lie 
algebras and of the algebra sl(2, R). These results enable construct- 
ing multi-parameter families of first- and second-order quasi-exactly 
solvable models. In particular, we have obtained two classes of quasi- 
exactly solvable matrix Schrodinger equations. 



1. Introduction 

There exists a small number of remarkable Hamiltonians (called exactly- 
solvable) whose spectra and corresponding eigenfunctions can be computed 

* e-mail : rzhdanov@ apmat . frecnct . kiev . ua 



1 



in a purely algebraic way (see, e.g. [|TJ). However, the choice of such Hamilto- 
nians is too restricted to meet numerous requirements coming from different 
fields of modern quantum physics. Recently, an intermediate class of Hamil- 
tonians was introduced by Turbiner and Ushveridze [[3] which allow an 
algebraic characterization of the part of their spectra. They call spectral 
problems of this kind quasi-exactly solvable. 

Quasi-exactly solvable models have an amazingly wide range of applica- 
tions in different fields of theoretical physics including conformal quantum 
field theories f|, solid-state physics || and Gaudin algebras (an excellent 
survey on this subject and an extensive list of references can be found in 
the monograph by Ushveridze 0). So it was only natural that there ap- 
peared different approaches to constructing quasi-exactly solvable models, 
including the one based on their conditional symmetries (for more details, 
see However, for the purposes of this paper the most appropriate is the 
Lie-algebraic approach suggested by Shifman and Turbiner p, [J and further 
developed by Gonzalez-Lopez, Kamran and Olver fl(|-|jl2|]. That is why we 



will give its brief description (further details can be found in ||). 

The Lie-algebraic approach to constructing quasi-exactly solvable one- 
dimensional stationary Schrodinger equations 

(-^ + V{x?Jil>{x) = \<il>(x) (1) 

relies heavily upon the properties of representations of the algebra sl(2, R) 

[Q , Q±] = ±Q±, [Q_, Q+] = 2Q 

by first-order differential operators. Namely, the approach in question utilizes 
the fact that the representation space of the algebra sl(2, R) having the basis 
elements 

Q_ = —, Q = x- -, Q+ = x- nx, (2) 

where n is an arbitrary natural number, has an (n+ l)-dimensional invariant 
subspace. Its basis is formed by the polynomials in x of the order not higher 
than n. Due to this fact, any bilinear combination of the operators (0) with 
constant coefficients yields a quasi-exactly solvable Hamiltonian 7i such that 
the equation TCip = Xtp can always be reduced to the form (|l|) with the help 
of a transformation 

^(x)^ F(x)^(f(x)). 
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Note that the above described procedure does not guarantee that eigen- 
functions of thus constructed quasi-exactly solvable Hamiltonians will be 
square-integrable. What can be done within this approach is to reduce a 
'differential' eigenvalue problem to a matrix one. The matter of a square 
integrability as well as other analytical properties of solutions obtained are 
to be investigated separately by independent methods (see, e.g. 

Recently, a number of papers devoted to constructing matrix quasi-exact- 
ly solvable models has been published |T^|-[pJ|. These papers use the same 



basic idea which is to fix a concrete subspace of sufficiently smooth multi- 
component functions and then to classify all second-order matrix differential 
operators leaving this subspace invariant. Furthermore, the above subspace 
is chosen to be the space of all multi-component functions with polynomial 
components. Posed in this way, the problem of constructing matrix quasi- 
exactly solvable includes as a subproblem a one of classifying realizations 
of Lie superalgebras by differential operators. Being very rich in interesting 
and important results this approach, however, contains an evident restriction 
which does not allow constructing all possible quasi-exactly solvable models. 
What we mean is the fact that an invariant subspace to be found is not 
necessarily formed by functions having polynomial coefficients. It is one of 
the results of the present paper that there exist principally different invariant 
subspaces. Thus there is a necessity for developing alternative approaches 
to the problem in question that do not require fixing a priori an invariant 
subspace. 

Our initial motivation for studying matrix quasi-exactly solvable prob- 
lems was to extend the list of exactly solvable Dirac equations of an electron 



via separation of variables. To this end we have suggested in |18| a method for 
constructing matrix quasi-exactly solvable models based on a direct general- 
ization of the Lie-algebraic approach for the case of multi-component wave 
functions. However, it turns out that the above method is universal enough 
to be applied for obtaining second-order quasi-exactly solvable models as 
well, including the Schrodinger equations with matrix potentials. 

Following we extend the class to which should belong basis elements 
of a Lie algebra under study (say, of the algebra s/(2,R)). We define this 
class as the set of matrix differential operators 

Q = Z{x)-£ + v(x), (3) 
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where £(x) is a smooth real-valued function, rj(x) is a smooth complex-valued 
r x r matrix function, and denote it as M. . The class M. is closed with respect 
to the binary operation 



and, consequently, form the infinite-dimensional Lie algebra. 

We will classify inequivalent representations of low dimension [d < 3) Lie 
algebras by operators belonging to M.. Next, we will study the additional 
constraints on the form of basis operators of the one- and two-dimensional 
Lie algebras imposed by the requirement that their representation spaces 
contain finite-dimensional invariant subspaces. These results will be used 
to obtain an exhaustive description of inequivalent representations of the 
algebra sl(2, R) by matrix differential operators (^) with r = 2. Composing 
linear and bilinear combinations of basis elements of sl(2, R) with constant 
matrix coefficients will yield multi-parameter first- and second-order quasi- 
exactly solvable matrix models. 

2. Classification of representations of low di- 
mension Lie algebras 

Since our aim is to get a quasi-exactly solvable model, we have to impose an 
additional restriction on the choice of the basis elements of the Lie algebras 
to be considered below. Namely, it is supposed that there exists at least 
one basis element such that the coefficient of does not vanish identically. 
This constraint is required to avoid purely matrix representations which are 
useless in context of quasi-exactly solvable models. 

Consider a first-order differential operator Q = £(x)d x + i](x) with £ ^ 0. 
Note that hereafter we denote ^ as d x . Let the function f(x) be defined by 
the relation 



and the matrix r x r function F(x) be a solution of system of ordinary 
differential equations 



{Qi,Q2}-^QiQ 2 -Q 2 Qi 



dcf 



[Qu Q2] 




fa) 



dF(x) 



+ rj(x)F(x) = 



dx 
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with det F(x) 7^ 0. Then the equivalence transformation 



Q^Q = (F(x))- 1 QF(x) 



with a subsequent change of the dependent variable 

/(*) 



x 



reduce the operator Q to become Q = d x . Consequently, any one-dimensional 
Lie algebra of first-order matrix differential operators Q = ^(x)d x +rj(x) with 
£ ^ is equivalent to the algebra (d x ). 

Abstract Lie algebras of dimension up to five have been classified by 
Mubarakzyanov in [19]. Below we give the lists of (non-zero) commutation 
relations which determine inequivalent Lie algebras of the dimension up to 
three. Note that the algebras which are direct sums of lower dimension Lie 
algebras are skipped from the lists. 



L 2> 1 


[Qu 


Q2] 


= Qu 










£3,1 


[Q 2 , 


Qs] 


= Qu 










£3,2 


[Qu 


Qs] 


= Qu 


[Q2, 


Qs] = 


Qi + 


Q2; 


£3,3 


[Qu 


Qs] 


= Qu 


[Q2, 


Qs] = 


Q2; 




£3,4 


[Qu 


Qs] 


= Qu 


[Q2, 


Qs] = 


-Q 2 




£3,5 


[Qu 


Qs] 


= Qu 


[Q2, 


Qs] = 


(1Q2, 


(0 < \a\ < 1); 


£3,6 


[Qu 


Qs] 


= ~Q2, 


[Q 


2, Qs] 


= Qx 




L3J 


[Qu 


Qs] 


= aQi - 


-Q2, 


[Q2, 


Qs]- 


= Q 1 + aQ 2 , (a > 0) 


£3,8 


[Qu 


Q2] 


= Qu 


[Qu 


Qs] = 


2Q 2 , 


[Q2, Qs] = Qs; 


£3,9 


[Qu 


Q2] 


= Qs, 


[Q2, 


Qs] = 


Qu 


[Qs, Qi] = Q2. 



Here a is a real parameter, the symbol L n ^ m stands for a Lie algebra of the 
dimension n numbered by m. 

Thus there exists only one two-dimensional Lie algebra L 2) i = {Q11Q2) 
which is not a direct sum of one-dimensional Lie algebras. 

If in operator Q\ = £,(x)d x + r](x) the coefficient £ is not identically 
zero, then using equivalence transformations defined at the beginning of this 
section we can reduce it to the form Q 1 = d x . Inserting Q 1 = d x , Q 2 = 
^{x)dx + fj(x) into the commutation relation [Qi, Q 2 ] = Q\ and equating 
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the coefficients of the powers of the operator d x yield system of ordinary 
differential equations for fj(x) 




1 




Hence we obtain £ = x + Ci, fj — A, where C\ G R is an arbitrary constant 
and A is an arbitrary constant r x r matrix. Shifting when necessary the 
variable x by a constant Ci we may put C\ — and thus get Q 2 = + A- 
If the operator Q\ has the form rj(x), then by convention the coefficient 
of d x of the operator Q 2 does not vanish identically. Consequently, there 
exists an equivalence transformation reducing the latter to the form Q 2 — d x . 
Substituting Q 1 = fj(x), Q 2 = d x into the commutation relation of the 
algebra L 2; i and equating the coefficients of the powers of the operator d x we 
obtain the following equations for fj(x): 



fj(x) = Ae x . 

Here A is an arbitrary r x r constant matrix. 

Summing up, we conclude that the two realizations of the algebra L 2 ,i 



exhaust the set of all possible inequivalent representations of the algebra in 
question within the class of matrix differential operators A4. 

In a similar way we have obtained complete lists of inequivalent represen- 
tations of the three- and four- dimensional Lie algebras within the class M. 
which are given below. 



dfj 
dx 



whence 



(1) Qi = Ae"*, Q 2 = d x - 

(2) Q 1 = d x , Q 2 = xd x + A 



(4) 
(5) 



L 



'3,1 • 



Q 1 = A, Q 2 = a 
[A, B] = 0; 



Q3 = ed x + Ax + B, 



L 



'3,2 • 



Qi = Ae~ x , C 
[A, B] = -eA; 



Q 2 = ee~ x d x + (B - Ax)e 



Q 3 = d : 
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L 3 , 3 : Qi = e- X (ed x + A), Q 2 = e~ x (ad x + B), Q 3 = d x , 
[A, B] — eB — aA; 

£3,4 : 

(1) Qi = Ae~ x , Q 2 =e x (ed x + B), Q 3 = d x , 
[A, B] = -eA; 

(2) Qi = 4, Q 2 = A, Q 3 = xd x + B, 
[A B] = -A; 

£3,5 : 

(1) Qx = Ae~ x , Q 2 =e- ax (ed x + B), Q 3 = d x , 
[A, B] = -eA; 

(2) Q 1 = e - X (d x + A), Q 2 = Be~ ax , Q 3 = d x , 
[A, B] = eaB; 

L 3 fi : Qi = A cos a; + B sin x, Q 2 = Bcosx — Asinx, Q 3 = d x , 
[A, B] = 0; 

L 3 j : Qi = e~ ax (Acosx + Bsinx), Q 2 = e~ ax (B cos x - Asinx), 
Q3 = d x , 
[A, B] = 0; 

^3,8 : Qi = d x , Q 2 = xd x + A, Q 3 = x 2 d x + 2Ax + B 
[A, B] = B; 

L 3i9 : No representations. 

In the above formulae a is an arbitrary constant, e = 0, 1, and A, B are r x r 
constant matrices. 



3. Quasi-exactly solvable matrix models 

As a second step of an implementation of the Lie-algebraic approach to con- 
structing matrix quasi-exactly solvable models we have to pick out from the 
whole set of realizations of Lie algebras listed in the previous section those 
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having finite-dimensional invariant subspaces. 

Consider first the one-dimensional Lie algebra (d x ). A space with basis 
vectors fi(x), . . . , f n (x) is invariant with respect to the action of the operator 
d x if there exist complex constants Ay such that 



dfi(x) 
dx 



x 



for all i = 1, . . . ,n. Solving this system of ordinary differential equations 
yields the following expressions for unknown vector-functions f,: 



=1 k=l 



where A is the constant nxn matrix having the entries A^; Cjy are arbitrary 
complex constants; the symbol {A)^ stands for the (i, j)th entry of the matrix 
A and e l5 . . . , e r are constant vectors forming an orthonormal basis of the 
space R r . 

It follows from the general theory of matrices that the above formulae 
can be represented in the following equivalent form (see, e.g. [20]): 



f,(x) = E^ X E^" ml (^)e fc . (6) 

3=1 k=l 

Here aii, . . . , a m are arbitrary complex numbers with |aj| < |a i+ i|, the symbol 
V\-^ m \x) stands for an (n — m)th degree polynomial in x, 1 < m < n, 
k — 1, . . . , n. 

As each realization of the low dimension Lie algebras obtained in Sec. II 
contains the operator d x , their finite-dimensional invariant subspaces are nec- 
essarily of the form (|6|). In what follows we will obtain complete description 
of finite-dimensional invariant subspaces of the representation spaces of the 
representations of the two-dimensional Lie algebra L 21 given in (Q), (|5|). 

First we turn to the case (f|). Let us study the restrictions on the choice of 
the basis vector-functions (](]) imposed by a requirement that the correspond- 
ing vector space V n is invariant with respect to the action of the operator 
Q2 = Ae~ x . By assumption, there exist complex constants Sij such that the 
relations 

n 

Ae x fi = ^2 
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hold with i = 1, . . . , n. Hence we get 

m r m r 

e ^- l)x e vr 3 r\x)Ae k = e e pJSrW, 

3=1 k=l 3=1 k=l 



where 



1 ijk - 2^ Ijk 
1=1 



T,-Plik m] (x)Ae k = (7) 



Comparing the coefficients of e aiX yield that = a* + 1, i — 1, . . . ,m — 1 
and furthermore 

r 

fc=l 

for all i = 1, . . . , n. 

Let us choose the new basis of the space R r in such a way that the 
first s basis elements ei , . . . , e s are eigenvectors of the matrix A with zero 
eigenvalues, namely 

Aei = 0, % — 1, . . . , s. 

Given this choice of the basis, it follows from (0) that V\i^ (x) = 0, k = 
s + 1, . . . , r. Hence, we conclude that the remaining basis vectors e s+1 , . . . , e r 
satisfy the relations 

s 

Aei = E a u e i' i = s + l,...,r 

3=1 

with some constant a^. 

Thus, the most general n-dimensional vector space V n invariant with re- 
spect to the Lie algebra (d x , Ae~ x ) is spanned by the vectors 

m r s 

Ux) = e« E e U ~ l)x E ^k m \^ k + e- E ^{x)^ 

3=2 k=l j=l 

where a is an arbitrary complex constant, Vj^ (x), (x) are arbitrary 

(n — N)th order polynomials in x, i — 1, . . . , n. And what is more the matrix 
A is of the following form: 
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where A is an arbitrary constant s x (r — s) matrix. 

Now we turn to the representation (|J). It is necessary to investigate 
the restrictions on the choice of the basis vector-functions (|6|) imposed by a 
requirement that the corresponding vector space V n is invariant with respect 
to the action of the operator Q 2 = xd x + A. By assumption, there exist 
complex constants Sy such that the relations 

(*± + A)t> = ±SA 

hold with % = l,...,n. Inserting the expressions (^) into these equations 
and comparing the coefficients of e ajX x k we arrive at the conclusion that 
<X\ — ■ ■ ■ — a m = 0. With this restriction the formulae give the most 
general finite-dimensional invariant subspace of the representation space of 
the Lie algebra (d x , xd x + A) 



V ik 1] (x)e k , i = l,... 



n. 



k=l 

A detailed investigation of finite-dimensional invariant subspaces admit- 
ted by the three- and four-dimensional Lie algebras is in progress now and 
will be the topic of our future publications. 

In what follows we will construct examples of quasi-exactly solvable two- 
component matrix models based on representations of the Lie algebra L^^ — 
sl(2, R). The construction procedure rely upon the assertion below which is 
given without proof. 

Theorem 1 The representation space of the algebra sl(2, R) having the basis 
elements 

Qi = dx, Q 2 = xd x + A, Q 2 = x 2 d x + 2xA + B, (8) 

where A, B are constant 2x2 matrices satisfying the relation [A, B\ = B, 
contains a finite dimensional invariant subspace iff the matrices A, B are of 
the form 

1) A=(-S _° ? ), B=(° J), or (9) 

2) A=f-t »), B=( » ?V (10, 
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Here n, m are arbitrary natural numbers with n> m. 

Using the fact that the algebra in question has the Casimir operator 
C = Bd x + A - A 2 it is not difficult to become convinced of the fact that 
representations of the form (g), (|9[) are the direct sums of two irreducible 
representations realized on the representation spaces 

Tlx = (§i,2rgi, . . . ,£ n e\), TZ 2 = (e 2 ,xe 2 , . . . ,x m e 2 ), 

where e x = (1,0) T , e 2 = (0, 1) T 

Representations (|J), (|10|) are also the direct sums of two irreducible rep- 
resentations realized on the representation spaces 

IZi = (nei, . . . , nx J ei + jr 7 e 2 , . . . , nx n e\ + nx n ~ e 2 ), 
TZ 2 = (e 2 ,xe 2 , . . . ,x n_2 e 2 ). 

According to the scheme given in the Introduction to get a quasi-exactly 
solvable model we have to compose a linear combination of basis elements of 
a Lie algebra of differential operators whose representation space has finite 
dimensional invariant subspace V n . And what is more, coefficients of this 
linear combination are constant matrices of the corresponding dimension 
whose action leaves the space V n invariant. 

Consider first the representation (|8]), (§) with n = m. According to 
the above its representation space contains 2(n + l)-dimensional invariant 
subspace V 2n+2 spanned by the vectors , e 2 x J , j — 0, ...,n. A direct 
verification shows that V 2n+2 is invariant with respect to action of any 2x2 
matrix. Composing a linear combination of (H), @ under n = m with 
coefficients being arbitrary 2x2 matrices yields the following quasi-exactly 
solvable model: 

(A x + A 2 x + A 3 x 2 ) du ^ + (At - nxA 3 )u = Au. (11) 

ax 

Here A\, A 2 , A%, A4 are arbitrary constant 2x2 matrices, u(x) is a two- 
component vector-function. 

Provided n — m + 1, the representation (H), @ gives rise to a quasi- 
exactly solvable matrix model of the form 

(A + B x x + ^x 2 )^^ + (S 3 + B 2 ^4)u = Au, (12) 
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where A\ is an arbitrary 2x2 matrix and Bi,B2,B 2 are arbitrary upper 
triangular 2x2 matrices. 

At last, if n > m+1, then the representation (|8|), (g) yields a quasi-exactly 
solvable matrix model of the form (|12~ ), where both A\ and Bi,B 2 ,B 3 are 
arbitrary upper triangular 2x2 matrices. 

A similar analysis shows that the representation (|T0|) gives rise to a 
quasi-exactly solvable model 

{B 1 + B 2 x + B 3 x 2 )^^- + {B A + B 3 {2xA + B))u = Au, (13) 

where 

B*=( X * b : ), i = l,...,4 



Aj, fej being arbitrary complex constants. 

Needless to say that any linear matrix model obtained from one of the 
above quasi-exactly solvable models by a change of variables is in its turn 
quasi-exactly solvable. In other words, equations (|Tl])-(|T3|) are representa- 
tives of equivalence classes of quasi-exactly solvable models. Other represen- 
tatives are obtained via transformation of variables 

x — > x = f(x), u — > F(x)u, (14) 

where f(x) is a smooth function and F(x) is an arbitrary invertible 2x2 
matrix whose entries are smooth functions of x. 

In the same way second-order quasi-exactly solvable matrix models are 
constructed. In particular, taking a bilinear combination of the operators (||) 
with A, B being given either by (|) or by flloD 



3 3 

h y2 a ijQiQj + X! a jQji 

where aty, 014 are arbitrary real constants, yields two families of quasi-exactly 
solvable matrix models of the form 

Tiw = Au. 

By a suitable transformation ( ]T4] ) the latter can be transformed to become 
a matrix Schrodinger equation 

— + V(x) u = Au. 



dx 
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Here V{x) is the 2x2 matrix potential, whose explicit form depends essen- 
tially on the parameters a^, ai and on the integers n, m. Let us stress that 
the matter of hermiticity of thus obtained matrix potential is by no means 
clear and needs special investigation. 

With a particular choice of the matrices A, B 

_° f ), B- 

the well-known nine-parameter family of scalar quasi-exactly solvable Schro- 
dinger equations is obtained [@-||. 




4. Conclusions 

The present paper is aimed primarily at solving the problem of classification 
of quasi-exactly matrix models by purely algebraic means. As a result, we 
get some classes of systems of first- and second-order ordinary differential 
equations such that a problem of finding their particular solutions reduces to 
solving matrix eigenvalue problem. Now to decide whether a given specific 
matrix model is solvable within the framework of an approach expounded 
above, one has to check whether it is possible to reduce it with the help of a 



transformation (|14|) to one of the canonical forms given in Sec. III. When one 
deals with a scalar model this check is being done trivially (see, for details, 
M). However for the case of matrix models it involves tedious and cumber- 
some calculations and is by itself rather nontrivial algebraic problem. As an 
illustration we will adduce an instructive example. Consider the following 
two-component matrix model: 

Hu = (iboiQi + iaa 2 Q2 + Cia% + c 2 a 2 ) u = Au, (15) 

where a,b,ci,c 2 are arbitrary real parameters with ab ^ 0, and <7i,a 2 are 
2x2 Pauli matrices, and 




Qi = d x , Q 2 = xd x + 

This model is quasi-exactly solvable by construction. Making the change 
of variables 

x = -sinh(ay), 
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w(x) = (cosh(a?/)) 1 / 2 exp I (c\ arctansinh(a?/) + C2 lncosh(ay)) 

x exp{— ia 3 arctansinh(m/)} ip(y) 



we reduce (15) to the Dirac-type equation 



iaM + o 2 V{y)J{y) = \^ 1 (16) 
ay 

where 

a 2 C2 — b 2 ci sinh(ay) 
abcosh(ay) 

is the well-known hyperbolic Poschl- Teller potential. This means that the 
Dirac equation (|1^) with the hyperbolic Poschl- Teller potential is quasi- 
exactly solvable. 

Thus an application of the obtained results to decide whether a given 
model is quasi-exactly solvable requires a considerable experience in manip- 
ulating matrix exponents. Generically, to check which equations of the form 



( |16l) can be reduced to one of the quasi-exactly solvable models constructed 
at the end of the previous section one has to solve systems of nonlinear 
algebraic equations. 

A technique used in the present paper can be generalized in order to 
enable one to classify multi-dimensional matrix models in a way as it was 
done for a scalar case by Gonzalez-Lopez, Kamran and Olver [TIJ-|r2| . 



The last important remark is that the property of quasi-exact solvability 
is intimately connected to the conditional symmetry of a model under study. 
This fact was for the first time noticed in our paper J/J, where we proved 
that quasi-exact solvability of stationary Schrodinger equations is in one-to- 
one correspondence with their conditional symmetry. As we believe, similar 
results can be obtained for matrix quasi-exactly solvable models as well. 
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